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Abstract. 

Physical processes influencing the properties of galaxies can be traced by the dependence and 
evolution of galaxy properties on their environment. A detailed understanding of this dependence 
can only be gained through comparison of observations with models, with an appropriate quan- 
tification of the rich parameter space describing the environment of the galaxy. We present a new, 
multiscale parameterization of galaxy environment which retains an observationally motivated sim- 
plicity whilst utilizing the information present on different scales. We examine how the distribution 
of galaxy (u-r) colours in the Sloan Digital Sky Survey (SDSS), parameterized using a double gaus- 
sian (red plus blue peak) fit, depends upon multiscale density. This allows us to probe the detailed 
dependence of galaxy properties on environment in a way which is independent of the halo model. 
Nonetheless, cross-correlation with the group catalogue constructed by Yang et al, 2007 shows that 
galaxy properties trace environment on different scales in a way which mimics that expected within 
the halo model. This provides independent support for the existence of virialized haloes, and im- 
portant additional clues to the role played by environment in the evolution of the galaxy population. 
This work is described in full by Wilman et al, 2010, MNRAS, in press. 
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INTRODUCTION 

The evolution of galaxies is heavily intertwined with the growth of the dark matter dom- 
inated structure in which they live, as baryons react to their local gravitational potential. 
This leads to a strong and evolving dependence of galaxy properties on their environment 
[e.g. 7, 3]. Environment is often defined in a observationally practical fashion, by count- 
ing the number of nearby neighbours. A different strategy is to assign galaxies to groups 
within a theoretical halo framework. Motivated by the need for a fair comparison of data 
with models, it would be preferable to define a statistical framework which retains the 
information contained by environment measured on different scales, without resorting to 
the inevitable assumptions required by a halo (group) finder during membership assign- 
ment. Simple two point galaxy statistics (such as correlation functions) trace the mean 
overdensity on different scales. However these are themselves closely correlated with 
one another. A better approach would be to examine how galaxy properties correlate 
with overdensity on one scale for a fixed overdensity on another, independently com- 
puted scale. Now that such a large, homogeneous sample of galaxies is available with 
the SDSS, this is now statistically possible [8, 4, 5, 6]. In Wilman et al, 2010, we take 




FIGURE 1. (a) Illustrates the annuli within which density is computed on different scales, as applied to 
three galaxies in different environments. On the three scales measured (starting with the smallest), these 
galaxies have 0,0,5 neighbours (lowest galaxy), 2,0,8 neighbours (middle galaxy) and 3,5,3 neighbours 
(top galaxy), (b) Demonstrating the bimodal distribution of SDSS galaxy (u-r) "model" colours. The 
colour distribution is well fit by a double gaussian distribution, described by five physically meaningful 
parameters (fred, ji red , Mbiue. Ored, Obiue)- 

the logical next step: defining a framework for quantifiable comparisons with models. 
The basic principles and initial results are described here. 

DATA AND METHOD 

Our approach is simple and provides a method to study the density-dependence of galaxy 
properties on different scales. For each galaxy the density of M r < —20. neighbours 
E ri ro is computed within a cylindrical annulus of depth dv rest ± lOOOkms 1 and in- 
ner,outer radii r ; ,r . This is computed on different scales, which, as illustrated in figure 1 
(a), traces different aspects of environment. 

Our sample is volume-limited within SDSS DR5: 0.015 < z < 0.08, M r < -20., pho- 
tometric completeness within 3Mpc of the galaxy > 98.645%, providing 73662 galax- 
ies. For each galaxy, on each scale, a spectroscopic completeness factor is computed 
(the fraction of potential spectroscopic targets with redshifts) and used to correct the 
measurements of density. 

To see how galaxy properties depend upon multiscale density, we turn to the bimodal 
(u-r) colour distribution of SDSS galaxies. It has been shown that this is well described 
by a double gaussian fit, and that the resultant parameters are strong functions of 
galaxy luminosity or mass and local density [2, 3, 1]. Figure 1 (b) shows an example 
of the double gaussian fit, illustrating the five independent parameters describing the 
distribution: the fraction of red galaxies f re d, the mean colour of each peak /i re( j, /ibiue 
and the widths of those peaks a re d, Obiue- Operating with the simplifying assumption 
that the red galaxy population is dominated by galaxies in which star formation has 
been truncated provides these parameters with relatively simple physical interpretations 



in terms of the overall star formation history for the galaxy population. 

RESULTS AND DISCUSSION 

Our sample is divided into bins, spaced evenly in log-density within a grid of two- 
scale density. These scales are always chosen so that the smaller scale measures the 
density of neighbours up to a fixed distance (dividing radius r^v), and the larger scale 
then measures the density in some annulus beyond r^y Figure 2 shows the fraction 
of red galaxies resulting from bimodal fitting within each bin for r div = 0.5,1.0 and 
2.0Mpc, where the outer radius for the outer annulus is 1.0,2.0 and 3.0Mpc respectively. 
There is also a (weak) dependence of f re d on luminosity within this sample, and so we 
repeat this procedure for each 0.5mag bin of luminosity in the range — 21.5 < M r < 
—20.0, with > 25 galaxies per luminosity and density/density bin required to produce 
reasonable fits (reduced ~ 1). For each bin we compute a residual density dependent 

A(f re d) = fred(M r , 

^0,r div ,£r div ,r o ) ~~ fred (Mr)- The three luminosity bins (which have 
similar appearance) are then coadded ( A(f ret j) is weighted by the number of contributing 
galaxies) to produce the final maps (figure 2). Colour gradients and contours illustrate 
the direction of change in A(f red ) (vertical contours illustrate no large scale dependence). 

We find a changing dependence of A(f re d) on density moving from small to large 
scales. On < 1 Mpc scales, there are more red galaxies in denser environments (a strong 
Eo ; o.5 dependence is supplemented with a residual positive correlation with Eo.5.i)- On 
larger scales there is no residual positive correlation, supporting previous results [8, 5, 6]. 
However on still larger scales > 2 Mpc, there is a highly significant anti-correlation with 
density at fixed < 2 Mpc scale density (i.e. more blue galaxies at higher £2,3)! 

The SDSS group (halo) catalogue of Yang et al. [9] is projected onto our multiscale 
density space. For the fraction of galaxies in haloes more massive than some fixed limit 
(e.g. 1O 12 - 5 M ), we find a dependence on multiscale density remarkably similar to 
that seen in figure 2. i.e. positive correlations on small scales (at larger density, more 
galaxies in massive haloes) are replaced by anti-correlations on large scales (at larger 
density, fewer galaxies in massive haloes)! More detailed analysis reveals the reason: For 
galaxies in the centre of haloes, density will always decrease with scale r, on average. 
Galaxies located a distance rh from the nearest halo centre will instead show increasing 
density with r up to r = r h (where the annulus encompasses the halo core), and then 
decrease to larger r. i.e. The galaxies at highest density at scale r, at fixed smaller scale 
density, have r l < r ~ r h < r . The virial radius of a 1O 12 5 M halo is ~ 0.2Mpc, so 
larger scales can trace the cores of haloes to which our galaxy does not belong. 

On larger scales this effect only becomes more important, leading to a lower fraction 
of galaxies in massive haloes at high densities on large scales - i.e. the anti-correlations. 
Therefore the higher fraction of blue galaxies at higher large scale density can be 
explained within the context of the halo model: These galaxies are in the accretion 
regions, and remain blue until after they have been accreted onto more massive haloes. 

This work is explained in more detail by Wilman et al., 2010, MNRAS, in press. Our 
multiscale density parameterization of the environment will provide an observationally 
motivated parameter space to examine the environment and the associated dependence 
of galaxy properties, and to make quantitative comparisons with expectations coming 




FIGURE 2. Dependence of the relative fraction of red galaxies, A(f re( j), on multiscale density. The 
bimodal fitting process is repeated independently for each bin in each two dimensional grid of interior 
and annulus density. This is illustrated for three choices of scales, moving from smaller to larger scales 
going from the top left to the bottom right. Grid colour changes and overplotted contours illustrate the 
direction of change within this parameter space. Interestingly our results illustrate a change from positive 
correlation with density (more red galaxies in denser environments) on small (< 1 Mpc) scales, to an anti- 
correlation with density (less red galaxies in denser environments) on large (> 2 Mpc) scales! This mimics 
the behaviour expected within a simple halo model framework (see text). 



from different physical models of galaxy evolution. 
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